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1. INTRODUCTION 

In order to predict wind turbine response charac- 
teristics in the presence of atmospheric turbu- 
lence, two major modeling steps are required. 

First, the important atmospheric sources for the 
force excitations TeTt by the wind turbine sys- 
tem must be identified and characterized. Second, 
a dynamic model must be developed which describes 
how thesia excitations are transmitted through the 
structure and power train. The goal of this 
paper is to establish the first modeling step, 
that of quantifying the important excitations due 
to the atmospheric turbulence. The dynamic model- 
ing cf the second step is undertaken in the ac- 
companying paper (1). 

Fluctuations in the aerodynamic forces on a wind 
turbine blade are generated by the relative motions 
of the air with respect to the blade. These rela- 
tive motions are comprised of two parts: the 

motions of the blade and the motions of the air. 

The motions of the air can further be divided 
Into the undisturbed turbulent flow and the 
"Induced. flow" due to the presence of t he w ind 
turbine vvake. The terms comprising the undisturbed 
flow will be characterized In this paper. More 
precisely, for a horizontal axis wind turbine, 
the aerodynamic forces are determined by the 
Instantaneous air velocity distribution along 
each of the turbine blades. These blades in turn 
are rotating through the turbulence field which 
Is being convected past the turbine rotor disc. 

It is thus necessary to characterize the wind 
turbulence field by a three-dimensional velocity 
vector which varies randomly with time and with 
the position In space. A complete statistical 
description of this turbulent velocity field re- 
quires the determination of all possible joint 
probability distributions between different 
velocity components at different times and 
positions in space. Clearly, such a description 
will not be possible without considerable simpli- 
fication. The validity of the resulting simpli- 
fied model will depend upon a comparison of the 
characteristics predicted by the model and those 
observed In the atmosphere and more importantly, 
those observed in actual wind turbine field tests. 
In this paper we will describe the assumptions 
and the analytical steps used to arrive at the 
simplified model. In the accompanying paper the 
model is used to predict wind turbine response 
characterti sties. It is hoped that these results 

will be verified In the near future by direct 
comparison with the results of actual field tests. 

2. MODEL ASSUMPTIONS AND APPROXIMATIONS 

The first assumption relates to the type of sta- 
tistical information which is necessary to describe 
the net aerodynamic forces and moments acting on 
the turbine rotor. Several authors (2,3) have 
indicated that the quantities needed for wind 


turbine design can be obtained from the mean and 
second-moment statistical characteristics of the 
various system responses. For stationary pro- 
cesses this information Is contained in the mean 
and power spectral density. In this type of 
analysis, the mean and power spectral density are 
characterized by a set of parameters. Rice's 
theory (4) for computing the frequency of level 
crossings or peaks is then used with the observed 
parameter probability densities to obtain the 
desired response statistics. In this paper, we 
will strive to determine the power spectral 
density characteristics of the turbulence. When 
they are combined with the machine dynamic model, 
we will assume that the resulting response sta- 
tistics will be useful for machine design. 

The next simplification assumes that the variation 
in the turbulent velocity observed at a stationary 
point is due primarily to the convection of the 
tu rbul ent eddies past the tower. Known as Taylor's 
frozen field hypothesis (5), this assumption is 
widely used in reducing fixed-tower, wind turbu- 
lence data and correlating these results with 
data from spatially separated points (6). 

The following assumptions which are often used in 
analyses involving aircraft flying through turbu- 
lence are more questionable when applied to turbu- 
lence observed in the atmospheric boundary layer. 
First, when the mean velocity field is subtracted 
from the total instantaneous velocity field, the 
resulting turbulent velocity is assumed to be 
locally homogeneous. Thus, when vertical separa- 
tions between points are as large as the disc 
diameter, the correlations are not explicitly 
heigiht dependent. Second, the field is assumed 
to be isotropic for all separations for which it 
is homogeneous. The latter assumption is known 
not to be precisely correct since the variance of 
the vertical component is less than the horizontal 
components (7) and the vertical and downwind com- 
ponents are correlated due to the boundary layer 
shear of the mean flow (8). However, no model 
currently exists for predicting the three-dimen- 
sional, nonisotropic correlations between velocity 
components at points separated in space. In the 
absence of a better model, the isotropic model will 
be used with the understanding that the results 
may need adjustment when more complete experimental 
results are available. 

With the previous assumptions {and assuming in- 
compressible flow), Batchelor (9) has shown that 
the correlation tensor between velocity components 
at spatially separated points has the form 

R^jd) =a2[f(c)6.j +icf(S)(5.j -^)] (2.1) 

where R^j(l) = E[v^. (x + |)Vj(x)] 

v^'(x) = ith velocity component at position'x 
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= ith component of the separation f 



f( 5 ) = longitudinal correlation function 

2 

o = variance of the turbulent velocity 
components 

Von Karman (10) suggested the form for the longi- 
tudinal correlation function 

f(C) ■ b(L)'« K, „(£,-) {2.2) 

where a = 1.339 

b = 0.5925 

ao 

L = integral scale - / f(c)dC 

0 

Ki/i(') = modified Bessel function of order 
1/3 

This function results in Kolmogorov's (11) -5/3 
power law for the inertial subrange in the longi- 
tudinal power spectral density. 

At this point, a very useful approximation due to 
Etkin (12) is introduced. The power of this 
approximation is that it separates the computation 
of the aerodynamic responses into two tractable 
pieces. In the first, the spatial variation of 
the turbulence is locally approximated by an ex- 
pansion. The various time varying turbulence com- 
ponents are then multiplied by standard aerody- 
namic influence coefficients to obtain the re- 
quired aerodynamic responses. These influence 
coefficients are the same as those that would be 
computed in the absence of turbulence. The re- 
sults of this procedure are extensively used in 
aircraft response calculations for flight through 
turbulence (13). The results for the airplane 
case, however, cannot be applied directly to the 
wind turbine problem because of major differences 
in the geometry. The aerodynamic surfaces of an 
airplane lie in a nearly horizontal plane while 
the blades of a horizontal axis wind turbine lie 
in a vertical plane nearly perpendicular to the 
mean wind. It is necessary then, to rederive 
the results in a form which is compatible with 
the wind turbine geometry. 

3. DERIVATION OF THE TURBULENCE MODEL 

The coordinate definitions used in this paper are 
shown in Figure 1. In the vadnity of the rotor 
disc, the turbulent velocity is expressed local- 
ly by the approximation 

vyCr,e,t) = Vyt) + v^y(t) rsinO + V^^^(t) rcos 8 

+ higher order terms (3.1) 

In this approximation, the spatial randomness of 
the turbulence is accounted for by the time vary- 
ing random quantities V-j(t), V^^x(t), 
higher order terms. While thi $ ’approximation ap- 
pears to be a Taylor series expansion, it is not. 
Because of the random nature of the spatial vari- 
ations, the samples from the statistical ensemble 


do not have the usual continuity and differentia- 
bility properties necessary for a true Taylor ex- 
pansion. The expansion, however, can be thought 
of as a functional approximation. Here, the ob- 
ject is to choose the terms in the expansion so 
as to minimize some measure of the approximation 
error. When dealing with random functions, a 
reasonable error measure is its variance. It will 
be understood that convergence of the approxima- 
tion series means that the error variance ap- 
proaches zero as more and more terms are included. 
Convergence in variance further implies that the 
series converges in probability (14), i.e., 

Lim Pr{|e |>£) = 0 for all e > 0 

n-K» 

where e is the approximation error including only 
the nth^order terms. 

At any given time, the terms V-j , Vi,x» and Vi,z 
are chosen to minimize the criterion 

€=jf - v.)^ dA (3.2) 

” A i ’ ’ 

where A = rotor disc area 

Vi “ ^i ^i X z 

= V^(r,e.t) 

The necessary conditions for the minimization are 
If- =f/U. -v.)dA = 0 

= f / {^ 1 - - V.) rsina dA = 0 (3.3) 

W . — "" f ^ dA = 0 

which in turn require that 

v.,(t) = 1 / v.(r,e.t)dA 

'^1 = V/ v.{r.e,t) rsine dA (3.4) 

A ’ 

V,- _(t) = f Vj (r,9,t) rcose dA 

i,z A ’ 

2 

where A = ttR the disc area 

I = I 7 = the area moments about the 
^ X and z axes. 

Thus, if the statistics of the turbulence field 
are known, then the statistics of Vj , V-j and 
Vi^z and any higher order terms can be determined. 
For example the autocorrelation function for the 
uniform, through- the-di sc component is expressed 
as 

R^^(t) ^ E [Vy(t+T)Vy(t)] 

= ^ / / E[v (r, 0 ,t+x)v (p,(f>,t)]dA,dA„ 

(3.5) 
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Using Taylor’s hypothesis yields 

Ry (t) * / Rop (Ci (3*6) 


curves can be plotted. These curves will be a 
one parameter family depending on the ratio of 
turbine size to turbulence integral scale (R/L). 
Example curves are shown in Figures 3 and 4. 


where C'l ” rsinO - psin{|) 

r = V T 

’2 \r 

= rcose - pcos4) 

dA-j = rdrd0 

dA^ = pdpd(|) 

and = mean wind speed 

In the isotropic case, 222 

= ohfU) + lcf'{c)(^^^)] (3.7) 

2 2 2 2 
where C = + ^3 

= - 2rpcos(e-(l)) + 

Even for the simple exponential correlation func- 
tion ^ 

f(d = e’ ^ (3.8) 


it is coubtful that an analytical expression for 
Ryjx) exists. Hence, numerical integration pro- 
cedures were employed to perform the required 
computations. Details of these procedures are 
found In the Appendix. 

At this point, it is convenient to rearrange the 
gradient terms for the in-plane components. This 
form is chosen because the resulting terms natu- 
rally appear when the velocity is expressed in 
components which rotate with the turbine blades. 
These terms can be Interpreted as local fluid ro- 
tations and strain rates. Thus, the following 
terms are defined 


Y 


xz 


i(V 

r z,x 




swi rl 


^xz * '^x.z^ 

^X2 2^^z,z ■ '^x.x^ 


shear strain rates 


(3.9) 




dilation 


Typical fluid streamlines giving rise to positive 
terms are shown in Figure 2. 


Retaining the uniform and gradient terms In the 
expansion results in* the following nine terms 
which vary randomly with time: Vj^, Vy, V^* Vy,x> 

Vy,z* Yxz. Yxz> Cxz» ^xz« The correlation sta- 
tistics of these terms can be computed using 
double-area integral expressions similar to 
Eq. (3.6). Because of the statistical isotropy, 
it is easily shown that all nine terms are 
mutually uncorrelated. Thus, all second moment 
statistics will be determineii by the autocorre- 
lation functions or the power spectral densities 
of the nine terms. Using the scaling parameters 
in Table 1, nondimensional power spectral density 


Table 1. Scaling Parameters 
Curves. 

for Nondimensional 


Seal ing 

Variables 

Parameter 

Turbulent velocity, V^. 

a 

Velocity gradient, V. . 

a/R 

Frequency, w 

V /L 
w 


Also shown In Figures 3 and” 4 "are approximate 
spectra derived from an exponential autocorrelation 
function. These approximate spectra match the 
computed spectra at low frequency and have the 
same total variance. Stationary, random processes 
with exponential autocorrelation functions can be 
conveniently represented by stochastic differen- 
tial equations of the form 


X + ax = bw 


(3.10) 


where x = random process 

w = white noise with flat PSD = q 

The autocorrelation function and power spectrum 
are 

Rx(t) = (3.11) 


S (ii)) — 2' 

a +u 


(3.12) 


respectively, from which the parameters a and b 
can be determined 


2RJ0) 

2R^(0) 


b = 




(3.13) 

(3.14) 


For the wind turbulence it is convenient to choose 
the white noise, power spectral density 

''IP 

W 

Nondimensional parameters can thus be defined 


(3.15) 


3 


(3.16) 
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R (0) 
2{-^) 


L b ^ 


uniform terms 


S (o)) = 2 2 
^ a + 0 ) 


4. MODEL ERROR DISCUSSION 


R^R (0) 

2(— I ) 

g 


(3.T7) 


gradient terms 


which win depend only on the ratio R/L. As an 
example of the computational procedure, consider 
the turbulence component Vy. In this case. 


1 E 

A— = -T / / g(f . 0) dA,dA„ 


(3.18) 


and V S (0) „ ” E ,'^w, 

^ ^ / / / g(f , -r-)dA,dA„(i^)dT 

La*^ A“^ 0 A A '■ 

„ r2 w2 2 

V T E “V T 

where g(^ , -^) = f(C) + |cr'(C)( ^ — ) 


^ = / r + - 2rpcos(e-())) + V^t 


dA^ = rdrdG 
dA^ = pdpd(|) 

and f(0 is the isotropic correlation function. 

The results of numerical computations for these 
Integrals are shown in Figures 5-8 for all of 
the turbulence components. 

In sunmiary, each of the turbulence terms are 
rm)deled by stochastic differential equations of 
the form 


X = ax + bw (3.20) 

where x = instantaneous value of one of the 

terms . ... . etc. 

w = nondimensional white noise with power 
spectral density q = a^L/V^ 


for uniform terms 


(3.21) 


(3.22) 


b* for gradient terms 

The nondimensional terms a* and b^ are found from 
Figures 5-8 as appropriate and depend on the ratio 
of turbine size to turbulence scale (R/L). Power 
spectral densities can be obtained if desired 
from the equation 


Three levels of approximation are introduced in 
this paper. In the first, the turbulence is 
modeled as locally homogeneous and Isotropic with 
correlations given by the Von Karman model. This 
assumption probably introduces the largest amount 
of error in the model. Several authors (15,16) 
indicate that the horizontal velocity components 
have a variance which is approximately three times 
the variance of the vertical component. If we 
assume that the turbulent velocity predicted by 
the isotropic model has a vertical component which 
is /3^ times too large, but is otherwise statis- 
tically correct, then the velocity error magnitude 
introduced has a variance 


e E [ E (5. - Vj^] 
i=1 ^ ’ 


= E [(V^v, - 




= 0.18 a*" (4.1) 

2 

where a = variance of horizontal components. 

The second level of approximation occurs in trun- 
cating the higher order terms in the spatial expan- 
sion, Thus, at any point on the rotor disc, the 
turbulent velocity is approximated by 


v^(r.0.t)=V.(t)+V^.^j^(t)rsine+V^^^(t)rcose (4.2) 


Since the velocity component through the rotor, 
Vy, produces the greatest aerodynamic force, con- 
sider the error variance produced by the approxi- 
mation of this component 


E^(r,e)= E[(Vy(r.0,t)-Vy(r,0,t))^] (4.3) 

Averaging over the rotor disc gives 


/ £.|(r,0) dA (4.4 

A 

Using the relations for terms Vy(t), Vy x(t) and 
y z(t) given by Equations 3.4 yields tfie useful 


/(v - V )v dA = 0 

y 

and hence 
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e, = -U / E[v^^]dA - -Jy f E[v V JdA 
1 Aa A ^ Aa'^ A ^ ^ 


Note> for a steady, uniform v^ocity, the force is 
constant and given by 


- 1 


-i ECVP ' 7T 


Aa" 


\ E[V 2](4.6) 
Aa*^ 


Now, let the approximate force be the result of 
the uniform and gradient terms in the model. Thus, 


which finally gives 

Ry (0) 

E] « 1 I ^ 

cr 


rV (0) 

Zi25 

2 



f(t) = ^ / r 


R'’ 0 


where = Vy(t) + 


Vy{r,! 2 t,t) dr ( 4 . 9 ) 

rsinnt + V (t)rcosJ 2 t 
y >2 


This quantity can be interpreted as a measure of 
the total variance of the part of the turbulent 
velocity that is not included in the model. Thus, 
tlii_iyjraged error, r, is zero if the approxima- 
te ^ _ 

and F = 1 if the trivial approximation v = 0 is 
usedj ^ 

In a similar fashion the quantity rg can be 
defined when only the uniform terms are retained 
and F2 when uniform, gradient, and quadratic terms 
are retained. Table 2 shows the effect of increas- 
ing rotor size relative to the turbulence scale. 


Integrating along the blade yields 
f(t)=C[Vy(t)+|| R(Vy^^(t)sinnt+Vy^^(t)cosnt)](4.10) 

The relative error variance Is given by 

^ - E[(f - f)^1 
’ E[f^] 

^ E [f^ ] - 2 E[ff] + E [ f^] (4 

E[f^] 


Table 2 - Relative 

Approximation 

Error Variance 

R 

L 

^0 

"1 

^2 

■ : _ ,01 

.044 

.026 

.023 

.054 (Mod M) 

.135 

.081 

.070 

f .1 

.201 

.121. 

.105 

,3 (Mod G) 

.397 

.250 

,218 

.5 

.527 

.348 

.304 

1.0 

.724 

.527 

.465 

2.0 

,889 

.737 

.663 


Observing the results given in this table, a sig- 
nificant Improvement is obtained when the gradient 
terms are included along with the uniform term. 
However, only a small improvement is obtained when 
the quadratic terms are also Included. This leads 
to the conclusion that the unmodeled portion of 
the turbulence Is highly disorganized and probably 
has a negligible effect on the forces and moments 
felt at the hub. 

To Investigate this effect further, the following 
aerodynamic model was assumed for a light, rigid 
blade cutting through the turbulent velocity field 


Substituting Equations 4.10 and 4.8 into these 
variance terms gives 


E[f^]=C^[Ry (0)+(^)^(Ry {0)sin^i^t+Ry (O)cos^nt)] 
y y.x 

( 4.12 

? ^ ^ — 

E[f ]=^ / / >"pqR2.r2jj,^2_p2j R22(Ci.O,53)drdp 


(4.13) 


where ^ = (r-p)slnJ^t 
= (r-p)cosf^t 

{•,-,•) = turbulent velocity correlation 
function 


E [ff] = C^pQ + 1^) 


(4.14) 


where 

^0 " n { ’" **22^^1 *0>E3)dAdr 


I, = -K f f r /p2 2 p cos (4i-nt)R,p(5T ,0,5.:,)dAdr 

I 4Rb 0 A R -r 22 1 3 


p 

f{t) = ^ / r /p2 ^2 v„(r,nt,t)dr (4.8) 

rJ 0 y 

where f = the net blade force (torque or 

thrust deviation from nominal) 

C = the aerodynamic influence coef- 
ficient 

u = rotation rate of the rotor 
and V (r,0,t) = instantaneous turbulent velocity. 


C-j = rsinQt - 

= rcosnt - pcoS(j5 
dA = pdpd(f> 

The normalized error variance, eg, defined by 
neglecting the gradient terms in the approximate 
velocity, is determined in a similar way. Table 3 
shows the results of these computations 
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Table 3 - Relative Blade Force Error Variance 


R 

L 

''o 

""l 

.01 

.024 

.009 

.054(Mod M) 

.076 

.030 

.1 

.116 

.046 

.3(Mod G) 

.248 

.104 

.5 

.350 

.157 

1 .0 

.536 

.278 

2.0 

.753 

.488 

Comparing the results 

of Table 

3 and Table 2, it 


Is seen that only half the unmodeled velocity vari- 
ance 1 s observed as unmodeled force variance. 

This result is due to the averaging effect of the 
Integration along the blade. If the blade were 
more realistically modeled with inertia it is ex- 
pected that little of the remaining unmodeled 
variance would be transmitted to the hub. 

The third level of approximation involves the use 
of the stochastic differential equation (Equation 
3,20) to model the uniform and gradient turbulence 
components. The accuracy of this approximation 
depends on how close the spectral form 


yields the stochastic differential equations 


+ b^W., 

(5.3) 

+ OXj + b2»2 

(5.4) 

82 X 3 + b 2 W 3 

(5.5) 


Since the original white noise inputs are uncor- 
related with identical power spectral densities, 
it can be shown that w^ and W 3 are also uncor- 
related white noise processes with the same power 
spectral density. This yields the following; matrix 
form for the stochastic differential equatiohs 


{X} = [A]{x> + [B]{w> 
{f} = [C]{x> 


where the matrices are given by 


[A] 

[B] 

[C] 


r-'i 

0 

0 ' 


0 

-"2 

Q 


0 

-n 

-^2~ 



0 

0 


0 

*^2 

0 


^ 0 

rr r 

0 

37T „ 

b2J 

m 



(5.6) 


(5.7) 


is to the spectra computed by Integration. Figures 
3 and 4 show two examples of such a comparison. 

The parameters a and b are chosen so that the total 
variance and the low frequency spectrum for the 
model are correct. 

Considering the results of these error calcula- 
tions, it is reasonable to expect that the turbu- 
lence inputs described statistically by the model 
will approximate the effect of the true turbulence 
on the wind turbine^ Realistic evaluation of the 
modeling error, however, can only be accomplished 
by comparison with experimental data. It is hoped 
that such a comparison can be made in the near 
future. 

5. AERODYNAMIC FORCE ON ROTATING WIND TURBINE BLADE 

As an illustration of how the turbulence interacts 
with a rotating turbine blade, consider the pre- 
vious example of a rigid blade rotating in the 
turbulent velocity field. Using the approximate 
turbulence model, the blade force is given by 


Using these equations the output power spectral 
density is given by the well known expression (17) 

S^(o)) = [H(iw)][Q][H^(-iui)] (5.8) 

where the row matrix of transfer functions is 
given by 

[H(iu))] = [C][iwl - A]’^[B] (5.9) 

Since the elements of the noise vector are uncor- 
related and have identical power spectral densities 


[Q] = 


q 

0 

^0 


0 

q 

0 


O"' 

0 

q ^ 


which gives 

S^(w) = q[H(iaj)][H'*^(-iu))] 

= q ^ 


(5.10) 


(5.11) 


f(t)=C[Vy(t)+^ R(Vy^j^{t)sinnt+Vy^^(t)cosJ5t)] (5.1) 


Defining the three components of the dynamic state 
vector 


x.|(t) = Vy(t) 

x„(t) = cosnt V (t) + sinot V (t) ) 

x^(t) = - sinfit V (t) + cosQt V^, ^(t) 

'5 J J 


(5.2) 


For the case at hand, 

Cb, C R b,(a„n-oo) 
[H(io>)] = [^. 2 2 

ai+ • 


(a2'^i(i3) 


C R 
(a2+iu))^+si^ 


and 


(5.12) 


(Cb-| ) q (^ ^ (dr,+Q +<i) )q 

■ -irr * (5.13) 


a-| + 0 ) 


(a2+f^ +0) ) -(2f^o] 
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Using the non-dimensional terms defined in Equa- 
tions 3.16 and 3.17 yields the result 


V^S^(<u) b*2)^(a*2 + 

" 7TTT 77TT~27 2 s 2 „ T2 

(3.14) 

where the non-dimensional frequencies are defined 
by 

0 >* = and £i* = ^ (3.15) 

The non-dimensional, power spectral density from 
Equation 3.14 is plotted in Figure 9, The para- 
meters for these blades were selected to corres- 
{)ond to two typical wind turbines of vastly dif- 
ferent size. Table 4 provides the key parameters 
for these two turbines. 

Parameters for Typical Wind Turbines 


Mod M Mod G 


Radius, R (ft) 

16.67 

150 

Rated Power 

8 kW 

2.5 MW 

Windspeed, (m.p.h.) 

16.63 

20 

Rotation Rate, Q (rpm) 

73.35 

17.5 

Turbulence Scale, L (ft) 

300 

500 


It is clear from Figure 9 that the effect of blade 
rotation is to concentrate the variance due to the 
turbulence gradient components at a frequency 
equal the rotation rate. This effect can be under- 
stood by considering the blade to be slicing 
through a slowly varying velocity gradient. As the 
blade encounters the higher velocity on one side 
of the rotor disc the force is increased. As it 
moves through 180'' the force reaches a minimum 
giving a fluctuating force at the rotor frequency. 
The importance of this effect can be seen by com- 
paring the relative contributions of the uniform 
and gradient components to the total variance of 
the blade force. Table 5 shows these results. 

Table 5 - Relative Contributions to 


Blade Force Variance 



Uniform 

Gradient 


Term 

Terms 

Mod M (8 kW) 

96% 

4% 

Mod G (2.5 MW) 

85% 

15% 


Clearly for the larger blade, 15 % of the variance 
at the relatively high rotor frequency could cause 
more fatigue damage than the 85% for the uniform 
component at the lower frequencies. 

6. CONCLUSIONS 

In this paper, we have formulated a theoretical 
model for the wind turbulence as it affects hori- 
zontal axis wind turbines. The model includes the 
effect of variations In the turbulent velocity 


across the rotor disc. An indication of the ap- 
proximation error in the model has also been given. 
It is expected that the model will be useful for 
determining how important the different turbulence 
effects are for given machine responses. This type 
of study has been made in the accompanying paper 
0). While we believe that the model wil l give 
qualitatively correct results, it Is important that 
experimental verification and any necessary model 
adjustments be made before it is used for design 
purposes. 
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APPENDIX 

Numerical Procedures Utilized in Model 
Development 

A,1 Area Integration Over Rotor Disc 

Two Gaussian quadrature formulas (18) were 
utilized to perform area integrations over 
the rotor disc. The distribution of points 
is shown in Figure 10 and given in Tables 6 
and 7. 


Table 6 - Sixteen Point Formula 


r . 
1 


.21132487 

.19634954 

.78867513 

.19634954 

Table 7 - Sixty-four 

Point Formula 

r . 

c. 

1 

1 

.26349923 

.03415057 

.57446451 

.06402420 

.81852949 

.06402420 

.96465961 

.03415057 


The quadrature formulas have the following 
form: 

n 4n 

/ f{r,e)dA = E z C.f{r.,0.) 

A 1=1 j=l ^ ^ ^ 



A six-point formula was also developed to re- 
duce the computational load. In this case, 
the radius was adjusted until the best match 
between computations using the sixty- four- 
point formula and the six-point formula was 
achieved and resulted in r = 0.69. For all 
computations, comparison was made between two 
formulas to verify accuracy to within five 
percent. 
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A. 2 Integration Along Radius 

An eight-point Gaussian quadrature formula 
(19) was utilized for radial integrations 
where required without angle dependence. 

1 n 

/ f(r)dr = E C.f(r.) 

0 i=l ’ ’ 

The weights and absissas are given in Table 


Table 8 - Radial Quadrature Formula 


'■l 


.0198550717 

,0506142681 

.1016667612 

.1111905172 

.2372337950 

.1568533229 

.4082826787 

.181 3418916 

.5917173212 

.181 3418916 

,7627662049 

.1568533229 

,8983332387 

.1111905172 

.9801449282 

.0506142681 


A. 3 Fourier Transforms 


Table 9 - Semi-Infinite Interval 


Quadrature 

Formula 




s 



.087649 

.206152 



.46270 

.331058 



1.14106 

.265796 



2,12928 

.136297 



3.43709 

.473289 

E-01 


5.07802 

.113000 

E-01 


7.07034 

.184910 

E-02 


9.43831 

.204272 

E-03 

* 

12.21422 

.148446 

E-04 


15.44153 

.682832 

E-06 

i 

19.18016 

.188100 

E-07 


23.51591 

.286235 

E-09 

- 

28.57873 

.212708 

E-11 

= 

34.58340 

.629797 

E-14 


41.94045 

.505047 

E-17 

- 

51.70116 

.416146 

E-21 



For calculation of power spectral densities 
the Fourier transform defined by 

S(u)) = / R(t)(Jt 


was numerically computed from the autocorre- 
lation function using the finite approxima- 
tion 

where ( 0 |^ = 

The fast Fourier transform techniques of the 
IMSL (20) library routine FFTRC were uti- 
lized. Several different time intervals were 
chosen to give overlapping spectra over the 
different frequency decades. 

A. 4 Semi-Infinite Time Interval 

To calculate the zero frequency power spectra, 
a Gaussian quadrature formula (21) was uti- 
lized. 


/ e'^f(x)dx = I C.f(x.) 

0 1=1 ^ ^ 

The absissas and weights for the sixteen 
points are given in Table 9. 


A. 5 Von Karman Correlation Function 

Central to all of the previous numerical inte- 
gration procedures is the need to compute the 
integrand. In this case the integrand will 
involve the evaluation of 

f(.) - 

and 

g(y) = J yf (y) 

where K, .»(•) = Modified Bessel function 
order of 1/3 

These functions are evaluated by the subrou- 
tine VK developed by Holley and Bryson (22). 
The method utilizes spline interpolation and 
asymptotic expansion giving a result accurate 
to six significant figures over a wide range 
of arguments. 
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Figure 1: Rotor Disc Coordinates 






Figure 2: In-Plane Velocity Gradient Terms 
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QUESTIONS AND ANSWERS 
W,E. Holley 


From; 

Q: 

A: 

From: 

Q: 

A: 

From: 

Q: 

A; 


L,P. Rowley 

How would you propose verifying your model? 

i-fitctid to compccvs th& stat^istvooly TesuXts pT&diotQd by the Tuodet W'lth vesuXts 
gavned from a planar anemometer array. This aomparison can be accomplished using 
results developed in system identification theory. We also would like to compare 
model predictions with wind turbine field data, 

T.E. Base 

Does your turbulence model satisfy the conservation equation, i.e., continuity? 

Yes, the Von Karman correlation function satisfies the constraint imposed by the 
continuity eguation. To the degree that the model approx'vmate s the Von Karman 
correlation function, it also satisfies continuity , 

L. Mi randy 

Is the only spatial effect in your model due to rsin 0, rcos 0 terms (gradients 
y , V, depend only on time) or do you have a spatial correlation like Dr. Sundar? 
r , X i / z 

The spatial variation effect includes only rstn0 and rcosB effects. Other more 
disorganized spatial variations are ignored. These higher order terms are the 
source of the error discussed in the paper. 
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